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The signature of bound state formation on the lattice is of particular interest in this talk. In the 
finite volume, where all states have discrete energies, it is rather hard to distinguish between a 
bound state and a scattering state if the bound state were close to threshold, i.e. like a "loosely 
bound state". To study bound states in the finite volume, we calculate the positronium spec- 
troscopy in the Higgs phase of U(\) gauge dynamics, where the photon is massive and then 
massive photons give rise to the short-ranged interparticle force. We try to identify bound state 
formation on the basis of the Liischer's finite size method, which suggests specific volume de- 
pendences of the energy gap/shift from the threshold energy for either bound states or scattering 
states. 
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1. Introduction 

Recently, a series of hadronic resonances have been discovered in various experiments [j|]. 
However, some of newly discovered states have unusual properties, which are not well understood 
from the viewpoint of the conventional quark-antiquark or three-quark states. Lattice QCD can 
potentially answer the question whether those states are really exotic hadron states since lattice 
QCD spectroscopy has been progressing with steadily increasing accuracy in the past several years. 

We are especially interested in some candidates of the molecular bound state: the A(1405) 
resonance as a KN bound state, the /o(980) and a (980) resonances as bound states of KK, the 
X(3&72) resonance as a weakly bound state of DD* and so on. In particular, such states except 
for the A(1405) are very close to threshold so that they would be "loosely bound states" like a 
deuteron. In the infinite volume, the loosely bound state is well denned since there is no continuum 
state below threshold. However, in a finite box on the lattice, all states have discrete energies. 
Even worse, the lowest level of elastic scattering states appears below threshold in the case if 
an interaction is attractive between two-particles [0]. Therefore, it is hard to clearly distinguish 
between the loosely bound state and the scattering state in this sense. 

In this paper, we present numerical studies of the bound state spectrum in the finite volume. 
As a pilot study of hadron molecular bound states, we explore the positronium spectroscopy in the 
Higgs phase of U(l) gauge dynamics, where the photon is massive and then massive photons give 
rise to the short-ranged interparticle force. In this model, we can control bound state formation 
in variation with the strength of the interparticle force. We then consider the application of the 
Liischer's finite size method which is relevant for elastic scattering of two particles with the 
finite range interaction, in order to identify the signature of bound state formation on the lattice in 
the finite volume. 



2. Liischer's finite size method 



Let us briefly review the Liischer's finite size method [gj. So far, several hadron scattering 
lengths have been successfully calculated by using this method. Especially, the I=2n—K channel, 
where the interaction is repulsive, have been intensively studied by one of authors [J3|] . 

It was shown by Liischer that the 5-wave scattering phase shift 8q is related to the energy shift 
in the total energy of two particles in the center of mass system in a finite box ^ : 
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where p and L are the relative momentum of two particles and the spacial extent, respectively. In a 
L 3 box with the periodic boundary, the generalized zeta function £¥oo(s,q) = J^nez 3 ( fl2 ~~ l) s 
is defined through analytic continuation in s. The asymptotic solution of Eq. (2.1) around q = 0, 
which corresponds to the energy shift of the lowest level of scattering states, is given by 
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with c\ = —2.837297 and C2 = 6.375183. jJ. denotes the reduced mass of two-particles as /I = 
mA -m^jimA The scattering length is defined through = lim^o tan 8 (p)/p. Eq. (2.2) 
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tells us that the lowest level of elastic scattering states appears below threshold on the lattice if 
an interaction is weakly attractive (oq > 0) between two particles. This point makes it difficult to 
discriminate between bound states and scattering states on the lattice. However, it is worth noting 



that the large L expansion formula up to 0(L ) in Eq. ( |2.2[ ) has no real solution of ciq for the case 



AE < — 2\ c ^l 2 w hile a o i s always (negative) real for AE > 0. A lower bound AE > — 2 \ Ci \hl 2 
may be crucial to identify the observed state below threshold with whether the lowest level of 
scattering states or a bound state. 

The question naturally arises as to how bound state formation is studied through the Luscher's 
formula since the quantum scattering theory implements bound state solutions. Indeed, another 
type of the asymptotic solution of Eq. ( |2.1[ ) around q = — oo, which was found by Seattle group [j|], 
represents a solution of bound states. Intuitively, non-vanishing negative energy gap AE in the 
infinite volume implies that a bound state is formed between two particles. This indicates that 
q = — oo in the limit of L — > oo is responsible for bound state formation. According to Ref. [g], an 
exponentially convergent expression of the the zeta function is given for negative q 
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where Y! ne2 ? means the summation without n 
(real y) as L — ► oo, Eq. (2J) leads to 
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(0,0,0). Suppose that p 1 approaches — y 2 < 



cot 8q(p) =i at q 
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in the infinite volume limit. This is certainly interpreted as bound state formation because the S- 
matrix S = e 2lS °^ = ^ot ^(pj— / nas a P°^ e at P 2 = ~Y Z - Therefore, for the bound state, one can 



derive the large L expansion formula around q = — oo from Eq. (2.1 ) 
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An L-independent term — ^ corresponds to the binding 



where (/?cot5o 

energy in the infinite volume limit. We can learn from Eq. ( [2.5D that "loosely bound state" is 
supposed to receive the larger finite volume correction than that of "tightly bound state" since the 
expansion parameter is scaled by y, which is associated with the binding energy. 



3. Compact Scalar QED 

To explore the signature of bound state formation on the lattice, we consider a bound state 
(positronium) between an electron and a positron in the compact QED with scalar matter: 

Ssqed[£/,<J>,^] = SmPM + £ %D w [U] x ,/¥ y (3.1) 

sites 

which is the compact U(i) gauge theory coupled to both scalar matter (Higgs) fields <I> and fermion 
(electron) fields *P. The action of "U (1) gauge + Higgs" part is described by the compact U (1)- 
Higgs model: 

S AU [U,<t>] = -p £ 9t{i^-A£9t{*X M **+/i}, W 

plaq. link 
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Figure 1: Schematic phase diagram of the compact U (l)-Higgs model in the fixed modulus case. A star 
mark represents our simulation point as (j3,/i) = (2.0, 0.6). 

where the constraint \<& x \ = 1 is imposed. In this study, we treat the fermion fields in the quenched 
approximation. We also consider the ^-charged fermion by replacing U{\) link fields as 

U x , il ^U c l ll =Yl c l =l U x ^ (3.3) 

in the Wilson-Dirac operator D\y- 

Fig.[TJ shows the schematic phase diagram of the compact £/(l)-Higgs model. There are three 
phases: the confinement phase, the Coulomb phase and the Higgs phase. The open symbols 
and filled symbols represent the second-order phase transition points (E: the end point {/3,/z} = 
{0.8485(8), 0.5260(9)} [Q] and S: the 4-dim XY model phase transition) and the first-order phase 
transition points (T: the tricritical point ~ {1,0.36} and C: the pure compact U{\) phase 

transition j3 c ~ 1). The lines ET and TC represent the first order line. The line TS corresponds to the 
Coulomb-Higgs transition, of which the order is somewhat controversial in the literature because 
of the large finite size effect. In this study, we have fixed j8 = 2.0 and h = 0.6 (the Higgs phase) 
where massive photons give rise to the short-ranged interparticle force. In the tree level, the vac- 
uum expectation value of the Higgs field and the photon mass are interpreted as (0hi ggs ) ~ a 
and Mphoton ~ a~ i \/h/P respectively. 

4. Numerical results 

We perform numerical simulations for positronium spectra (^o and 3 S\ states) in the Higgs 
phase of t/(l) gauge dynamics on L 3 x 24 lattices with several spatial sizes, L = 8, 10, 12, 16, 
20, 24. Two-point functions of So and 3 Si states are constructed from the bilinear pseudo-scalar 
operator ^75^ and vector operator ^VxJ^x respectively. To evaluate a threshold energy, we also 
calculate the electron mass in the Landau gauge. 
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Figure 2: Masses of 3 Si and 'So positronium states as functions of the electron mass. The left (right) panel 
corresponds to results with charge-three (four) electrons. The dotted lines represent the threshold energy, 
which is evaluated by 2 x aM e i ec tron- 

Figs. ^|show masses of ! So and 3 Si positronium states as functions of the electron mass. The 
certain energy gap from the threshold energy appears in simulations with charge-four electron 
fields. It is natural to expect that higher charged electrons provide the larger energy gap since 
the interparticle force is proportional to (charge q) 2 . For q = 4, the hyperfine mass splitting of 
positronium is also clearly observed. 

The volume dependences of energy gaps at aMeiectron — 0.5 are shown in Figs. ||[ All data 
points in the right panel (q = 4) are clearly below the lower bound for the asymptotic solution of 
the scattering state. The volume dependence is drastically changed around L ~ 12 — 16. Data for 
the larger lattice sizes are well fitted by the form inspired by the asymptotic solution of the bound 
state, Eq. (2.5). The energy gaps for either 1 Sq or 3 Si states should remain finite in the infinite 
volume limit. Therefore, bound states of electron-positron are certainly formed in simulations with 
charge-four electrons even in the Higgs phase, where interparticle forces are short-ranged. 

On the other hand, an upward tendency of the L dependence is observed as spatial size L 
increases in the right panel (q = 3). However, all data points are located near the lower bound 
for the asymptotic solution of the scattering state. We also remark that the L dependence seems 
to become opposite around L ~ 20 — 24. These observations suggest that the observed states are 
unlikely the lowest level of elastic scattering states, rather likely a "loosely bound state". However, 
to make firm conclusions on this, more detailed study and also the data for larger L are required. 



5. Summary 

We explored the signature of bound state formation in positronium spectra ('So and 3 Si states) 
in compact scalar QED model, where the short-range interaction between electron and positron can 
be realized as in the Higgs phase. In the case of highly charged electrons (q = 4), the energy levels 
of both 'So and 3 Si are found to be far below threshold, while electron-positron states with the lower 
charged electron (q < 3) appear slightly below, but close to threshold. For q = 4, we found a specific 
volume dependence of an energy gap between the total energy of electron-positron states and the 
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Figure 3: Mass gaps from the threshold energy as functions of spatial lattice size L in lattice units. Dashed 
lines represent a lower bound for the asymptotic solution of the scattering state, which is given in the text. 
Solid curves in the right panel (q = 4) are fits of the form AEl = A£oo + j-e~ bL , which is inspired by the 
asymptotic solution of the Liischer's formula for the bound state iQ. 

threshold energy, which is well described by the asymptotic solution of the Liischer's formula for a 
bound state. More detail analysis, which includes the sensitivity test of mass spectrum with respect 
to spatial boundary conditions and the examination by the volume dependence of spectral weights, 
is now under way to investigate the formation of "loosely bound state" in the finite volume. 
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